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Summary 

This  paper  discusses  the  space  harmonic  Fourier 
analysis  of  magnetic  fields  for  a  compensated  pulsed 
alternator  (compul sator)  or  an  active  rotary  flux  com¬ 
pressor  (ARFC),  which  has  been  applied  for  two  particu¬ 
lar  sets  of  boundary  conditions  appropriate  for  these 
machines  to  calculate  magnetic  field  distributions, 
inductances,  and  how  those  inductances  vary  with  rotor 
position.  The  results  of  the  analytic  model  are  com¬ 
pared  with  experimental  results.  This  paper  also  shows 
the  results  of  the  extensive  parametric  studies  using 
the  inductances  calculated  from  space  harmonic  distri¬ 
bution  (SHD)  code  to  explore  design  alternatives  to 
improve  the  performance  of  these  machines. 

Introduction 

The  operation  of  a  compensated  pulsed  alternator 
(compul sator)  or  an  active  rotary  flux  compressor 
(ARFC)  involves  a  number  of  phenomena,  the  most  impor¬ 
tant  of  which  is  the  delivery  of  output  current  through 
a  circuit  of  varying  inductance  that  provides: 

a.  compensation  to  reduce  inductive  impedance 
to  current  flow,  and 

b.  flux  compression  by  changing  inductance  to 
enhance  output  voltage. 

Because  the  operation  depends  intimately  on  the 
magnetic  fields  in  the  machine  and  how  they  vary  in 
time  and  space,  it  is  important  to  have  as  accurate 
a  magnetic  field  analysis  as  practical.  The  basic  pro¬ 
blem  thus  becomes  one  of  calculating  inductance  as  a 
function  of  rotor  position. 

Both  the  armature  winding  and  the  compensating 
winding  occupy  known  radial  and  azimuthal  regions  in 
the  air  gap. 

A  solution  of  Laplace's  equation  for  a  two-dimen¬ 
sional  problem  in  cylindrical  coordinates  with  radial 
and  azimuthal  variations  is  a  sinusoid;  hence  if  the 
excitation  and  boundary  conditions  can  be  expressed  as 
Fourier  series,  then  the  fields  in  the  machine  can  be 
found  as  space  Fourier  series. 

Fields  and  Stored  Magnetic  Energy 

The  machine  geometry  to  be  analyzed  is  shown  in 
Figure  1.  The  origin  (^  -  =)  =  o  corresponds  to  the 
center  line  of  the  belt  of  the  windings  that  has  a 
missing  conductor,  which  is  the  result  of  the  serpen¬ 
tine  winding  configuration.  Both  armature  and  compen¬ 
sating  windings  are  shown  as  occupying  certain  radial 
and  azimuthal  regions  defined  by  the  angles  and  radii. 
If,  at  a  particular  radius  in  one  of  the  windings,  we 
plot  the  current  density  as  a  function  of  angle,  we  get 
the  current  distribution  shown  in  Figure  2.  Thus, 
assuming  a  series  of  the  form 
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Mg.  1.  Geometry  of  the  machine  used  for 
space  harmonic  distribution 
inductance  calculation 


winding  current  sheets 

interior  of  the  machine  where  the  conductors  carry  cur¬ 
rent  axially  and  end  effects  can  be  neglected,  a  differ 
ential  element  of  the  current  density  can  be  treated  as 
an  axial  (z  -  directed)  current  sheet,  which  for  the 

th  ,  . 

n  component  is 

Si  =  \  Kn  cos  n  ^  “  “)  (2) 


J  =  J  +  5  J  cos  n  (f-«),  (1)  where  K  =  J  dp. 

o  n=1  n 

Field  Intensities 

the  coefficients  are  found  in  the  usual  way. 

Using  the  radial  variable  p  to  denote  a  location  The  next  step  is  to  find  the  field  produced  in  the 

within  a  current  distribution  and  considering  the 
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machine  by  the  current  sheet  of  Equation  2  for  each 
winding  and  then  integrate  over  both  winding  dis¬ 
tributions  to  find  the  total  field  strength.  Because 
this  is  a  linear  field  problem,  each  Fourier  component 
for  each  winding  can  be  considered  separately  and  then 
the  results  may  be  summed  later. 

To  proceed  with  the  analysis  we  consider  only  a 
current  sheet  described  by  equation  2  in  an  air  gap 
with  free-space  permeability  u  .  Except  within  the 

current  sheet,  which  is  infinitesimal  (dp  )  in  thick¬ 
ness,  and  possibly  at  the  boundaries,  there  are  no 
currents  flowing,  in  which  case 

V  x  H  =  0,  (3) 

which  allows  H  to  be  described  as 


H  =  -vU.  (4) 

Because 

B  =  H  (5) 

and 

V.B  =  0  (6) 

the  potential  U  must  be  a  solution  of  Laplace's 
equati on 

v2U  =  0  .  (7) 


In  the  interior  of  the  machine  axial  variations  of 
fields  can  be  ignored;  thus  Laplace's  equation  for 
cylindrical  coordinates  can  be  written  with 

-r|-  =  0  to  obtain 

1  3  /  \  1  32U 

V2U  =  7  W  fr  jy.  j+  r*  W  -  0.  (8) 

The  variables  can  then  be  separated  by  assuming  a 
product  solution 

U  =  R  (r)  0  (e)  (9) 


substituting  it  in  Equation  8  and  solving  the  two  re¬ 
sulting  differential  equations  to  obtain  the  general 
sol  uti  on 


U  = 
n 


+  D 


A  sin  ne  +  B  cos  ne 


(10) 


For  the  general  current  sheet  of  Equation  2, 
which  is  located  at  r=  p,  two  regions  are  identified: 


(The  field  is  calculated  at  a 
region  i :  <r  <p  (point  radially  inwards  from  the 

(current  sheet. 

(The  field  is  calculated  at  a 
region  o:  p  <r  <R$  (point  radially  outward  from  the 
(current  sheet. 


Subscripts  i  and  o  will  be  used  to  denote  the 
region,  with  respect  to  the  current  sheet,  in  which 
the  field  is  being  calculated. 

The  two  sets  of  boundary  conditions  for  a  compul- 
sator  or  an  ARFC  are  shown  in  Table  I.  Both  boundary 
conditions  I  and  II  are  for  a  laminated  rotor  made  of 
ferromagnetic  material  of  high  permeability.  Boundary 
condition  I  pertains  to  solid  steel  poles  whereas 
boundary  condition  II  pertains  to  laminated  stator 
magnetic  material . 

The  application  of  these  two  sets  of  boundary 
conditions  and  the  resulting  field  expressions  are 
shown  in  detail  in  references  1  and  2. 


TABLE  I 

Two  Sets  of  Compulsator-ARFC  Boundary  Conditions 


At  Radius 

Boundary  | 

Condition  I  ; 

Boundary 

Condition  II 

Rotor  Surface 

r  =  \ 

Infinitely  permeable 

H  .  =  0 

ei 

Stator  Surface 

r  =  Rs 

Perfectly 

Conducting 

H  =0 
ro 

infinite  ly 
Conducting 

H  =0 

60 

At  Current 

Sheet 
r  =  p 

H  continuous 

r 

H  .  =  H 
r  i  ro 

At  Current 

Sheet 
r  =  p 

^Heno  "  Heni^ 

=  Kn  cos  n  (e  -  °=) 

To  calculate  the  total  field  due  to  single  winding 
consider  first  the  compensating  winding  which,  accord¬ 
ing  to  Figure  1,  occupies  the  radial  space 

Rc  <  p  <  Rc . 
a  b 


For  this  winding  three  regions  are  identified  in 
which  the  field  must  be  calculated  separately  as 
follows:  n 


for  Rr  <r  <Rg  Hnc  =  / d  Hnic 

Kr 


(n : 


for  R5  <r  <R6  Hnc 


r  r\c 

/dHnoc*/‘ 


d  "„ic 


(12) 


for  Rg  <r  <RS  H 


r 

■/< 


nc  Hnoc' 

Re 


(13) 


The  use  of  these  integrals  leads  to  the  field 
expressions  for  the  compensating  winding.  Because  the 
process  is  exactly  the  same  for  the  armature  winding, 
the  field  expressions  for  the  armature  winding  are  ob¬ 
tained  from  those  for  the  compensating  windings  by  mak¬ 
ing  the  replacements: 


cc 


c 


cc 

a 


References  1  and  2  contain  these  final  field 
expressions  for  the  two  sets  of  boundary  conditions. 


Magnetic  Stored  Energies 


The  next  step  in  this  space  harmonic  field 
analysis  is  to  use  the  field  intensities  to  calculate 
the  stored  magnetic  energy  and  from  that  the  self  and 
mutual  inductances  of  the  windings.  Basically  the 
magnetic  stored  energy  is  calculated  from 


/ 


1  pH2  d  vol . 

2  0 


(14) 


It  will  be  convenient  to  calculate  the  self-  and 
mutual  inductances  separately.  For  either  winding 
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(15) 


yields 


2ir 


w  = 

m 


/  /1//2  *J0(Hr  +  %)lr  dr  de- 


Because  different  harmonics  are  orthogonal  with 
respect  to  the  integration  on  6,  all  cross  terms  will 
integrate  to  zero  and  Equation  15  can  be  written  as 


M 


m 


2 

n=l 


M 

run 


(16) 


wi  th 

2T 

Wmn  -JJ  V2  uQ(H2n  +  H2n)fr  dr  de.  (17) 


Each  field  component  can  be  written  in  the  form 


fn  =  Hrnm  sin  n(e-a) 

(18) 

en  =  Henm  cos  n(9'a)- 

09) 

Squaring  these  expressions,  replacing  the  squares  of 
the  sinusoidal  functions  with  trigonometric  identities, 
and  carrying  out  the  integration  on  e  yields  zero  for 
the  sinusoidal  terms  leaving 

■nvnl  r 

W  =  — — -  /  (H2  +  H2  )r dr.  (20) 

mn  2  J  '  rnm  enirr 

Using  the  field  expressions,  the  stored  energies 
for  the  various  regions  for  armature  and  compensating 
winding  excited  separately  were  calculated.  The  re¬ 
sulting  expressions  are  tabulated  and  displayed  in 
references  1  and  2. 

Using  the  fact  that  for  self-inductance 

Wm  =  1/2  Li2  (21) 

and  the  expressions  for  energies,  along  with  the  rela¬ 
tions  for  the  maximum  compensating  and  armature  wind¬ 
ing  current  densities  yields  the  self-inductance  for 
use  in  the  expression 


L  = 


5 

n=l 


(22) 


The  final  calculation  to  be  made  is  of  the  mutual 
inductance  between  the  armature  and  compensating  wind¬ 
ings.  This  requires  the  calculation  of  magnetic  stor¬ 
ed  energy  resulting  only  from  simultaneous  excitation 
of  both  windings.  From  Equation  15  these  cross  terms 
yiel  d 

Wm  =  //  uo^HrcHra  +  HecHea^r  dr  d6‘  <23> 
o 


Using  the  orthogonality  of  different  Fourier  components 
and  the  general  form  of  the  field  expressions  given  in 
Equations  18  and  19  with  the  nth  component  along  with 
the  use  of  trigonometric  identities  on  the  products  of 
the  sinusoidal  functions  and  integrating  on  e  yields 


W 

mn 


■nvol 


i 


(H  H  +  H  H  )cosn(a  -a  )rdr. 

'  rnmc  rnma  enmc  enma'  '  a  c '  ^4) 


n  f  2 


M 


n 


5 

2 

k=l 


Wmnk 


n  = 


2 


5 

2 

k=l 


m2k 

1aic 


(27) 


(28) 


When  the  armature  and  compensating  windings  are 
connected  in  series,  the  total  inductance  of  the 
circuit  is 


L  =  L  +  L  +  2M  (29) 

d  C 

which  with  the  use  of  Equation  26  becomes 

L  =  La  +  Lr  +  2  2  COSn(ci  -a  )  (30) 

a  c  n_i  n  a  c 

which  describes  explicitly  how  the  circuit  inductance 
varies  with  the  relative  angular  position  (a  -a  )  of 
the  two  windings.  3  c 

Compulsator-ARFC  Inductance 

The  calculated  magnetic  field  strengths,  magnetic 
energy,  and  effective  inductance  for  the  two  sets  of 
boundary  conditions  were  used  to  develop  computer  pro¬ 
grams  called  space  harmonic  distribution  (SHD)  codes. 
Two  machines  that  match  our  first  set  of  boundary 
conditions,  viz,  laminated  rotor  and  solid  stator 
poles,  had  been  designed  and  fabricated  at  the  Center 
for  Electromechanics  at  the  University  of  Texas  at 
Austin  (CEM-UT).  These  were  the  engineering  proto- 

1  3 

type  compulsator  and  the  desk  model  compulsator  . 

The  second  set  of  boundary  conditions,  viz,  laminated 
rotor  and  laminated  stator,  was  an  attempt  to  increase 
the  magnetic  flux  compression  ratio  of  a  compulsator 
or  an  ARFC.  An  ARFC  was  designed  and  fabricated  with 

laminated  rotor  and  laminated  stator  at  CEM-UT.4  The 
optimum  angular  span  for  the  windings  and  conductor 
configuration  were  selected  with  the  help  of  the  SHD 

p 

code  corresponding  to  boundary  condition  II.  The 
exact  winding  configuration  and  system  geometry  of  the 
three  machines  were  readily  available  to  use  in  calcu¬ 
lating  the  effective  inductance  of  each  machine.  For 
boundary  condition  I,  since  the  poles  are  not  perfect 
conductors  there  will  be  some  penetration  of  the  field 
into  the  pole  region.  An  effective  Hr  =  0  boundary 

(assuming  y  =  yQ  in  the  region  between  Rs  and  Rj)  was 

located  at  a  certain  depth  into  the  pole  which  yielded 
SHD  code  results  that  matched  the  measured  values. 

The  calculated  and  measured  values  of  effective  mini¬ 
mum  and  maximum  inductance  of  the  three  machines  are 
given  in  Table  II . 

In  using  the  SHD  code  the  effect  of  the  end  turns 
on  the  calculated  inductance  per  unit  length  must  be 
considered  unless  the  machine  has  many  poles,  in 
which  case  the  length  of  the  end  turns  is  small  com¬ 
pared  to  the  active  length  of  the  rotor. 


The  use  of  the  field  expressions  with  this  equation 
results  in  the  mutual  energy  storage  expressions  in  the 
five  regions.  Combining  these  energy  expressions  with 
the  maximum  compensating  and  armature  winding  current 
density  expressions,  the  expression 


W _ =  i  i . 


mn 


n 


and 


M  =  2  M  cosn (a  -a  ) 
n=l  n  3  c 


(25) 

(26) 


Parametric  Design  Studies 

The  two  SHD  codes  for  the  two  sets  of  boundary 
conditions  were  used  to  study  the  variation  in  induc¬ 
tance  as  different  parameters  are  changed  in  the 
engineering  prototype  compulsator  and  the  ARFC.  Only 
the  parametric  design  study  for  the  prototype  compul¬ 
sator  will  be  discussed  here.  The  parametric  design 
study  for  the  ARFC  is  discussed  in  detail  in  reference 
2. 
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TABLE  II 

Calculated  and  Measured 
(LCR  bridge  measurement  at  1.0  kHz) 
values  of  effective  minimum  and  maximum  inductance  (H) 


'"max 

L  . 
min 

Flux  com¬ 
pression  ratio 

Machine 

calc 

meas  ! 

calc 

meas 

cal  c 

meas 

Proto¬ 

type 

compul¬ 

sator 

177. 7n 

1  76y 

24.1m 

26.3m 

7.37 

6.69 

Desk 

model 

compul¬ 

sator 

138u 

136u 

16.4m 

17.5m 

8.39 

7.77 

ARFC 

1.19m 

1.09m 

26.4m 

23.2m 

45.08 

46.98 

Self-Inductance  vs.  Radial  Position  of  the  Winding 

To  determine  the  sensitivity  of  the  self-induct¬ 
ance  of  the  armature  winding  to  the  inner  radius  of 
the  winding  the  SHD  code  was  suitably  modified  and  the 
calculations  for  the  prototype  compul sator  were  repeated 
over  the  range  Rr  <_  <_  R$  -  (Rg  -  R^ ) .  The  self¬ 

inductance  of  the  armature  winding  varied  from  80mH  to 
1 OyH  as  the  inner  radius  of  the  armature  winding  was 
increased  from  the  rotor  surface  until  the  outer 
radius  coincided  with  the  pole  faces. 

Effective  Inductance  vs.  Radial  Separation  of  Windings 

The  flux  compression  ratio  decreases  signifi¬ 
cantly  as  the  radial  separation  ( Rg-R^ )  is  increased 

To  achieve  a  narrower  pulse  width  it  is  necessary  to 
operate  the  device  at  the  smallest  air  gap  that 
mechanical  and  electrical  constraints  will  allow. 

Effective  Inductance  vs.  Circumferential  Span  of 
Windings 

The  circumferential  distribution  of  the  windings 
is  important  in  determining  the  effective  inductance 
as  a  function  of  the  angular  position  of  the  rotor 
winding  relative  to  the  compensating  winding  during 

the  pulse.  The  calculated  values  of  (Lm  /L  .  )  are 

max  min 

plotted  against  compensating  winding  conductor  span  at 
N=12  turns  (e^)  in  Figure  3.  The  radial  thickness  of 

the  windings  is  assumed  to  remain  constant. 

To  maintain  the  same  winding  resistance  the  same 
variation  of  conductor  span  is  considered  but  with  a 
constant  conductor  cross-section.  Flux  compression 
ratio  decreases  as  the  conductor  span  is  reduced  for 
the  constant  cross-section  case.  Thus  a  wide,  radially 
thin  conductor  design  is  preferable  for  maximum 
compression  ratio. 

Number  of  Poles 

The  pulse  width  of  the  discharge  can  be  decreased 
by  increasing  the  number  of  poles.  But  as  shown  in 
Figure  4  the  flux  compression  ratio  decreases  with 
increasing  number  of  poles. 

Conclusions 

The  detailed  space  harmonic  analysis  described  in 
this  paper  enabled  the  successful  establishment  of  ex¬ 
pressions  for  fields,  stored  energy  and  inductances 


Fig.  3.  Flux  compression  ratio  of  the  prototype  compul 
sator  versus  conductor  span--constant  radial 
thickness 


Fig.  4.  Flux  compression  ratio  vs.  number  of  poles  in 
a  1-m  diameter  laminated  rotor,  laminated 
stator  ARFC 

that  agree  well  with  measurements  from  three  experi¬ 
mental  models.  The  analytic  model  provides  a  versa¬ 
tile  tool  for  parametric-design  studies  of  both  a 
compulsator  and  an  ARFC.  The  corresponding  SHD  code 
facilitates  the  calculation  of  the  space  harmonic 
current  densities,  magnetic  stored  energies,  self- 
and  mutual  inductances,  and  the  flux  compression 
ratio.  Over  and  above  proving  the  viability  of  the 
compulsator  concept,  the  SHD  technique  for  calculating 
effective  inductances  is  a  useful  design  tool. 
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